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I. INTRODUCTION 



In the last four decades there has been an extraordinary progress in solving the Einstein- 
Maxwell equations with two Killing vectors The most important solutions have been 
related to charged black holes, i.e., Reissner-Nordstrom static black hole and Kerr-Newman 
rotating black hole. In addition, it is well known the solution of Plebanski-Demianski (PD) 
lU metric, containing the last solutions as particular cases; the PD solution is the most 
general solution in the Petrov classification type D. Besides, recently, these solutions were 
generalized in the metric affine gravity 0-^]. 

The present work is the result of an investigation that we did using computer algebra 
methods of the family of solutions of Einstein-Maxwell equations with cosmological constant 
A 

Gi,u + Xqi^u - KTf,y = 0. (1) 



In particular, we have extensively used REDUCE 3.6 



The Maxwell equations can be written in tensorial form as 

dfsF^^ = r, (2) 

d^Ff,e + deF^p + dpFe^ = 0, (3) 
where F"^ is the Maxwell tensor 

Fab = dbAa - daAb, (4) 
with the quadripotential A°', and the current density quadrivector (in this work we take 

r = 0). 

II. STRUCTURE OF THE PD METRIC 

Following reference we consider the class of space-times with signature (— , -|-, -|-, -|-) 
in real generalized coordinates Q^P^ cr) with the line element given by 

ds' = gooix^dr' + 2go3{xndrda + guix^dq^ + ^?22(x")rfp' + gssixlda^ (5) 

where we take only the following two-coordinate dependence of = X^'iQ^p)- The PD 
metric has a similar structure [Q 

ds' = ^ - P'dcrr + ^dq' + ^dp' + '^idr + q'da)'] , (6) 

or 

ds = -jj^dr + ^ drda + j^dq + ^dp + da . (7) 

The structure functions V and Q depend on the coordinates p and q, respectively. The 
functions H and A are defined as H = 1 — pq and A = p^ + q^, respectively. This space-time 
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has two Killing vectors 8/ Or and 8/ da that commute between them. The vector d/dr is 
timelike in the region where Q — V > 0, and the space-time is stationary there. 
The real electromagnetic potential is given by 



A" = ^ ((eoq + gop), 0, 0, (gop - eop)pq) , 



(8) 



where eo and go are the electric- and magnetic-like charges, respectively. 

It is well known that the structure functions V and Q are polynomials of fourth degree, 
whose coefficients, depending on seven parameters, are classified as mass m, NUT parameter 
n, angular momentum jo, acceleration parameter b, cosmological constant A, and the charges 
eo,go- 



V = (h-gl)+ 2np - ep2 + 2mp^ - b + e^ 



Q = (b + eg) - 2mq + eq^ - 2nq^ 



3^ 



P 



b-g^ + -lq^ 



(9) 
(10) 



where 6 — 7 — g. 

We mention that recently the following generalization of this metric has been found by 
Garcia and Macias Q, including a new parameter /x (new acceleration parameter) 



V = 


(b - g^) + 


2np — ep^ + 2m/ip'^ 


Q = 


(b + e^) - 


2mq + eq^ — 2n/xq^ 


H = 


1 - /ipq, 




A = 


2 1 2 
p +q . 





/^'(b + eg) + ^ 



(11) 

(12) 

(13) 
(14) 



III. FAMILY OF SOLUTION FOR PD METRIC 

We present now the family of solutions for Einstein-Maxwell equations with cosmological 
constant A 

G^u + ^gfiu - /tT^i. = 0. (15) 

Using the program for Excalc package [|TI| in REDUCE 3.6, we can introduce the fol- 
lowing ansatz for the polynomials 'P(p) and Q(q) 

^(p) = Po + PiP + P2P^ + PsP^ + P4P'^, (16) 

Q(q) = qo + qiq + q2q^ + qsq^ + q4q^- (17) 

Via computer algebra we obtained the following relations 
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qo-Po = eo + go, 



P3 

qs 

P2 
P4 



-/^pi, 

-q2, 
A 

'3 ~ 



/^^qo, 



q4 = P4 + /i^(eo + go)- 



(18) 
(19) 
(20) 
(21) 

(22) 

(23) 



At this point we generalize the relation (|T8]) by introducing two new real parameters, a and 
j3, as follows 



qo - Po 



h+ia + l)el + Pgl - h + ael + ip - l)g, 



el + -' 



0- 



(24) 



Finally, the family of solutions given via the structure functions V and Q, can be written 
by taking into account the structure function for the PD metric 



^(P) 



Q(q) 



h + ael + {p- l)gl\ + 2np - ep^ + 2m/ip^ 

A" 



/i'(b + (a + l)e^ + /5g^) + 



P 



b + (a + l)e^ + pgl] - 2mq + eq^ - 2n/iq^ 

A' 



/x'(b + «e^+(/?-l)go') + 



q 



H = 1 - ^pq, 
A = p2 + q2. 



(25) 



(26) 

(27) 
(28) 

1, a = 



The PD metric corresponds to the following values in the new parameters: fx 
and P = 0. 

This family of solutions have the property that they remain solutions under the inversion 
transformation 



1 



Q 



(29) 



performed in both the structure functions and the electromagnetic potential [|l|. In addition, 
we can find the cylindrically symmetric counterpart of the metric (^) in similar form to 
Garcia et. al, 131- 



IV. ROOTS AND SINGULARITIES IN THE PD GENERALIZED 

When we consider the problem of obtaining the roots and singularities in this generalized 
metric, via calculating the invariants R, Rq/jR"^, Rapes R"^^^ , we found an intrinsic singu- 
larity for A = + = 0. Under the coordinate transformation (q — ^ r, p — — jo cos 6'), the 
function A = gives the corresponding ring Kerr singularity, where the new parameters a 
and (3 do not modify this property. 
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The roots in the polynomials V = (|25|) and Q = ( p6D represent singularities in the 
coordinates, and therefore they are removable via an adequate coordinate transformations. 
The singularities in the coordinates, show us some defficiencies in the use of these coordinate 
systems to describe the space-time under consideration. 

The roots in the polynomial V = do not have any physical interest, and we omit them. 
However, the real roots in the polynomial Q = 0, represent null hypersuperfaces of event 
horizon type with q = constant. 



The solution for Q = is obtained via conventional methods [|T2[ . Consider the following 
form 



b+ {a + l)el + I3gl - 2mq + eq^ - 2n^q^ 



^l\h + ael + {P-l)gl) + 



X 



or, 



q + a^q^ + 02?^ + flig + Qq = 0, 



q' = 0, (30) 



(31) 



where 



as 
ao 



2nfi 



e 



A ' 

3 



/x2(6 + ae2 + (/3-l)(72) + | 
2m 

/i2(6 + ae2 + (/3-l)^?2) + |' 
[6 + (a + l)el + pgl 



A ' 



^^{h + ael + {P-l)gl) + ^^ 
Using the Ferrari method for equations of fourth degree and the Cardano method for those 



of third degree |jT2|, we can get the following roots 



where 



A 



B 



C 



(32) 



A^_2^(tl)l/3 + (i2)l/3_I 



f + f -a2 + (ti)^/^ + (t2)^/^ 



1/2 



{f-A) 

h 1 , to = U — ti, 

2 V 4 27/ ' 
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u = bi — 



fo2 
"2 



V = bo- 



6162 2bl 
^ ^ 27' 



60 
bi 



qi I- J 
{mnij, + Aqoq^} , 62 = —, 



fi\b + {a + l)el + pg^) + 



A 



The polynomial Q (]30|) can be rewritten as 

q4 jq'' + asq^ + a2q^ + aiq + ao} = q4(q - qi)(q - qii)(q - qiii)(q - qiv), (33) 

where qi, qn, qm and qiv are the roots (0) in ascendent order of magnitude, satisfying the 
following relations 



IV 



51 ?i = -03 

i=i 

J2 lilj = ^2 = 



2nfi 



^2 (5 + «e2 + _ 1)^2) + A 

e 



i>j 



i>j>k 

T^iqi = ao 



^^(b + ael + {P-l)gi) + l 
2m 



-ai 



^,^{b + ael + {P-l)gl) + l 
b+{a + l)el + Pgl 



^,^{b + ael + {(5-l)gl) + \ 



(34) 
(35) 
(36) 
(37) 



where ai, i = 0, 1, 2, 3 are the coefficients of Q. 

Also, we can obtain red-shifted hypersuperfaces for this generalized solution, which are 



given by the equation [13 



5^00 



P-Q 



0. 



(38) 



The roots of (|38|) are analogue in form to those in equations (0), except that the coefficient 
ao is now a function of the coordinate p. 



V. REDUCED FAMILY OF SOLUTIONS FROM THE FAMILY OF SOLUTIONS 

TO PD METRIC 



Certain family of reduced solutions can be obtained from the family of solutions of the 
PD metric when we drop appropriately the new parameters and perform the coordinate 
transformations from (r, p, q, o") (other coordinates) as used in Garcia [Q]. We give these 
reduced generalized metrics in table I. All of them correspond to stationary axisymmetric 
space-times with Maxwell field matter source. 
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Metric 
electromagnetic potential 



Metric 


structure functions 


PD generalized 
(10 parameters) 


p =[b + ael + {l3- l)g'^] + 2np - ep^ + 2m/xp3 

- [ii\h + (a + l)e2 + + A] p4 
Q = [6 + (a + l)eg + (3gl\ - 2mq + eq^ - 2niiq^ 
-[^,^{b + ael + {(^-l)gi) + l] q' 
A = p'^ + q'^ 
H = 1 — ^pq 


PD standard 
= l,a = 0,/3 = 0) 


-p = {h- gl) + 2np - ep^ + 2mp-^ - (6 + + ^)p'^ 

Q = \h + el) - 2mq + eq^ - 2nq^ - {b - g^ + 
A = p2 + ^2 

H = 1 ~ pq 


PC generalized 
(/. = 0) 


r = [b + ael + {p- + 2np - ep^ - Ap4 

Q = [b+{a + l)el + pgl] - 2mq + eq^ - \q^ 
A = p2 + g2 


KN generalized 
= 0, n = 0, A 7^ 


= bo' + «eg + (/? - 1)^0^] - (1 - - Ip' 
Q = \ji + (« + + Pgi] - ^mq + (1 - %q' - \q^ 
A = p2 + g2 

= 1 


KN generalized 
= 0,n = 0, A = 0, 


V = [jI + axi + (/? - l)^?^] - p^ 

Q = bo + (" + l)eo + /^^o] - 2mg + q^ 
A = + g2 

= 1 


KN standard 
^ = 0,n = 0, A = 0, 
a = 0,13 = 1, 
^ ^ Jo, e ^ 1 


r = jl-p' 

Q = (jo + eo + ^?o) - 2mg + q^ 

A = p2 + g2 

= 1 



Table I 

Generalization and standard cases contained in the PD solutions 



VI. WHAT TO DO THE NEW PARAMETERS a AND /3? 

In order to study the features of the new parameters a and /3, we take the KN metric as 
an example. As it is known, the exterior gravitational field of a black hole bounded by an 
event horizon is characterized by three parameters: mass m, angular momentum jo, and the 
electric charge cq. 

When we take the Boyer-Lindquist coordinate transformation in the generalized PD 
metric, we obtain the generalized KN metric (with A = 0) through the transformation 
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T : {t ^t- jo0, q^r,p^ -jo cos 9, cr 



.} 



(39) 



where now the space-time have the following coordinate = (^5^5^5 0)) the explicit 
form of the metric is 



ds 



2 _ T^K N - QkN ^_^2 



aA 



KN 



A 



fo + r')VKN - fo sin' ^Qxtv} dtd(t> + ^dr' 



Q 



KN 



+ + -2X 1 Uo + ) ^KN - Jo sm OQkn \ 

rKN Jo^KN ^ ^ 

with the structure functions Vkn, Qkn and Akn given by 



Vkn 
Qkn 



jI + (a + l)el + (3gl 



— 2mr + 



Aii-Tv = r + Jo cos 61. 

The hypersurfaces of infinite shift would satisfy the condition 

'Pkn — Qkn 

9tt = T = 0, 

^kn 

impliying that there exist two such hypersurfaces at the radial distance 

1/2 



where the equation 



I 2 / 2 I 2 I -2 2 

m ± m — (Co + (^o + Jo ' 



Q<el + gl+3l<m\ 



(40) 



(41) 
(42) 
(43) 



(44) 



(45) 



(46) 



gives us the condition for non existence of imaginary values and of negative radius. 

Examining the quadratic invariants, we find an intrinsic ring singularity for Akn = 
that is similar to Kerr's singularity 

AKN = r^+JoCos^O = 0. (47) 

As is usual, the event horizon can be calculated from the equation Qkn = 0, leading to 

1/2 



m ± 



m^-((a+l)e^ + /?(?o+Jo) 



(48) 



We can conclude that the parameters a and f3 modify the radius of the event horizon of the 
generalized solution of KN metric. We investigate this modification in the following. 



A. Parametrization of the event horizon in the KN space time 



The ideal case under consideration corresponds to the equation (^Sf) , where there are no 
event horizon with imaginary or negative radius. Therefore, we have the condition 

0<{a + l)el + Pgl+j^<m^ (49) 
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or 



- {el + Jo') < ael + (3gl <rr? - (e^ + j^). 



65 a 



Pi(0,m2 
P2(0,gg) 




a + go /? < 



eoa + go/5 = -(eo+jc 



(50) 



Fig. 1.: Family of solutions to the equation (Q). 

The set of solutions shown in Figure 1 represents a family of straight lines of slope 
7 = — 1. If the values for a and (5 given at the points P2 or P3 are introduced in the metric 
(^), the standard KN solution and the corresponding event horizons are recovered. In fact, 
any value for a and (5 coinciding with a straight line crossing through the points and P3 
will lead to the standard KN metric. Then, the family of parallel straight lines given by 
the equation (^Dj) represents different KN space-times for each case, each with its own event 
horizon. 

Any straight line that crosses the solution range of a and /5, intersecting the family of 
parallel lines of slope 7, will lead to all spacetimes found in this range of solutions. Then 
we suggest that there is a simplification in the set of parameters if we take the straight line 
as the axis a (/? = 0) and introduce a new parameter v that reruns this set of solutions. We 
propose the following form of this straight line 



2 / 2 I ■2\ 



< z/ < 1. 



(51) 



Thus, the new structure functions are 



V 



KN 



/ 2 I 2 I ■2-\ 



Qk at = m z/ — 2mr + r , 
AxTV = + cos^ 9, 



(52) 

(53) 
(54) 



where z/ satisfies equation (PH). 

All black hole properties mentioned above are mantained with the inclusion of this new 
parameter and the important point is that this solution, it not violate no-hair theorem. 

The equation Qkn = for the structure function (l5^ ) determines the event horizons 
whose radiis are given by 



m 



< z/ < 1. 



(55) 
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Therefore, from the range of v we get the following values for the interior r_ and exterior 
r+ radiis 

< r_ < m, m < r+ < 2m. (56) 

One can see that the corresponding solutions are proportional to the geometric mass. The 
results given in equations -(|5^) represent the parametrization for the event horizon of 
KN type black hole, where each value of v gives a different space-time with the interior and 
exterior horizon. 



B. Quantitave analysis to the KN space-times and its event horizons 



It is easy to see that with v = {e^ + + jQ)/m^ we can recover the standard KN 
space-time. 

Equation (^) gives us two ideal cases 



and 



< eo + ^o +Jo < m\ 



2 I 2 I -2 2 

eo + 9o+Jo = ^ ■ 



(57) 



(58) 



The set of solutions for the KN spacetimes taking into account the parameter z/, are shown 
in figure 2. 



eo + go+Jo 



eo + go 




m 



1±(1 



■2m 



■m 







Fig. 2.: . Set of solutions for the KN space-times for the case < Cg + ffo + Jq < m? ■ The superposed 
straight hues r_ and r_|_ represent the radius of the interior and exterior event horizon, respectively. 

To analyse this figure, we can divide it in two regions 



I:Q<u<{el + gl+3l)/m' 
II:{el + gl+3l)/m'<u<l. 



(59) 



The first region I is not physical, because the event horizons cross the infinite redshift 
surface in different points of the rotation axis, and the corresponding angles of intersection 
for both surfaces violate the range for the function sin^x (0 < sin^a; < 1). 

The region II describes the physical parametrization, because it contains the family of 
KN space-times with their event horizons of radiis 
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rKN± = rn[l ± Vl - H- 



(60) 



The set of spacetimes corresponding to the value = 1, represents non-extreme black holes 
with only one event horizon aX, tkn — fn. 

VII. CONCLUSIONS 

Although, wc dcah with a difficult algebra we have been able to obtain a generalization 
of the PD metric by extensive use of the REDUCE 3.6 system. The latter could be used in 
any available computer resources. 

We introduced a PD metric structure with 10 parameters and an appropriate reduction of 
parameters in the generalized PD metric leads to obtain two new families of solutions to the 
Plebanski-Carter and Kerr-Newman metrics, containing the new parameters a and j3. We 
point out that the introduction of the parameters a and (3 does not modify the generalized 
solutions when the electromagnetic charge parameters Cq and Qq are interchanged, as it does 
in the structure functions and the electromagnetic potential considering all their posible 
combinations. 

Interpretation of the generalized KN metric as the exterior field of a collapsed star 
surrounded by an event horizon, black hole, has lead to a new equivalent generalized 

metric, but now simplified with only one dimensionless parameter called v. In addition, we 
have determined the way v modifies the radius of event horizon such that the corresponding 
solutions are proportional to the geometric mass. On the other hand, the redshift surfaces 
remain constant and independent of v. 

The minimum value of the parameter v corresponds to the space-time of a standard KN 
type black hole. Increments in the value of v increase the region between the redshift surfaces 
and the event horizons (the ergospheres) until the maximum value for v — 1. This maximum 
value corresponds to the space-time of a KN type black hole with only one event horizon 
which is not necessarely an extreme black hole. The space-time of a standard extreme KN 
type black hole can also be recovered. 

Dropping the electromagnetic charge parameters allows to get a generalized Kerr solution 
that corresponds to a vacuum stationary axisymmetric solution of Einstein equations. The 
analysis and interpretation as a black hole and parametrization of the event horizons is 
completely analogous to the generalized KN case. 
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